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ABSTRACT: 

General theorem describing a relation between diagonal of one-electron den- 
sity matrix and a certain class of many-electron ensembles of determinant states 
is proved. As a corollary to this theorem a constructive proof of sufficiency of 
Coleman's representability conditions is obtained. It is shown that there exist rig- 
orous schemes for construction of energy of many-electron system as functionals 
of one-electron density. 
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Introduction 

In density functional theory (DFT) approaches it is accepted that the 
electronic energy of many electron systems (at least for the ground state) 
can be presented as a functional of the first order density, or, in other words, 
of the diagonal of the first order density matrix. In this connection the 
following question seems to be pertinent: Are there rigorous schemes, that 
do not involve approximations and hypothesis of any kind, for construction 
of electronic energy as a functional of the first order density ? In present 
paper we make an attempt to answer this question. 

The second question closely related to the first one may be formulated as: 
Is it possible to consider electronic energy of many-electron system for some 
fixed basis as a function of occupation numbers ? Recently orbital occupancy 
(00) approach developing this idea for DFT type functionals was formulated 
[1-4]. In present paper we show how energy expressions involving diagonal el- 
ements of the first order density matrix may be rigorously constructed. Such 
energy expressions may be used only in non-gradient optimization schemes 
since these expressions are not differentiable in classic sense. 



Basic Definitions 

For a fixed basis set of n orthonormal spin-orbitals the corresponding 
finite-dimensional Fock space Tn is spanned by determinants \L) where L 
runs over all subsets of the spin-orbital index set N. Its p-electron sector 
J-n, p is spanned by determinants \L) with \L\ = p. Basis determinants will 
be labelled by subsets and all sign conventions connected with their repre- 
sentation as the Grassman product of ordered spin-orbitals will be included 
in the definition of the creation-annihilation operators. 

The set of g-electron density operators is defined as 

S N , q = {t q € F N , q ® r NA : ij = i, & > & Tr{t q ) = 1} (1) 
The diagonal mapping over F N ® J 7 ^ is 

d(t) = J2(L\t\L)e L , (2) 

LcN 

where t G JFjy ® J 7 ^ and e L = \L)(L\. 

The contraction operator over J-n®J~n is defined in terms of the standard 
fermion creation-annihilation operators as 

n 

c = <g> a\. (3) 

i=\ 

Definition 1. q-electron density operator t q G £jv,g is called weakly 
p-representable if there exists p-electron density operator t p G £n, p such that 

^(d(t p )) = d(t q ). (4) 

This definition is correct because the contraction operator possesses the prop- 
erty 

c{d(£ N , p )) C d(£ N , q ). 

The set d(£jv, P ) is called the standard (unit) simplex of the operator space 
d(^N, P <S> J~N, P ) an d its characterization is given by 

(p) (p) 
T N>P = d(£ N , p ) = {J2 ^e L : A L > & ^ A L = 1}. (5) 

LcN LcN 
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The combinatorial structure of 1n, p is very simple: Any part of the set 
of all p-element subsets of the index set N determines a face of T n p and its 
complementary part generates the opposite face. In particular, there are (™) 
hyperfaces opposite to the corresponding vertices. 

Definition 2. 

W N , p , q = ^ c p ~id(S N>p ). (6) 

Weak representability problem may be formulated as the problem of de- 
scription of the polyhedron W^, p , q with arbitrary admissible n,p, and q. 
Since, by definition, density operators are Hermitean, this polyhedron may 
be conveniently embedded into the real Euclidean space Mat, 9 of the dimen- 
sion ( n ) with its canonical basis vectors labelled by g-element subsets of 
N. With such an embedding the tensor products of the fermion creation- 
annihilation operators involved in the expression (3) should be replaced by 
the commuting (Bose) annihilation operators 

0jej \o, if j i j ' {{) 

acting on the vector space 

n 

R N = R Ntq . (8) 

WN,p,q is a polyhedron situated in the real afline hyperplane 

(?) 

H Ns = {\e R N , q :J2 X J = ^ ( 9 ) 

JCN 

Let us put 

! 1 (q) 

w plq (L) = q -^e L = -j- eK. (10) 

\q) KcL 

Directly from definition it follows that the polyhedron WN, P , q is the convex 
hull of (™) vectors w p [ q (L): 

W Nj>tq = Conv({w plq (L)} LcN ). (11) 

To the best of our knowledge, in contrast to the parametric description given 
by Eq.(ll), the analytic description (that is the description in terms of the 
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hyperfaces) of this polyhedron is obtained only for the case q — 1 and is given 
by the following assertion that is just a consequence of the general theorem 
by Coleman [5-10]: 

Theorem 1. Polyhedron Wn, p ,i is the set of solutions of the system 

< A fe < i k G N 

E A, = I (12) 

jeN 

This polyhedron has 2n hyperfaces with normals 

v° k =pe k , (13a) 

and 

vl = -pek + ^ej, (136) 
where k G N, and e k are canonical basis vectors of the Euclidean space W 1 . 

Restoration of p-Electron Wave Functions from One-Electron 
Density Matrix Diagonal 

With arbitrary vector A^ G W^, p ,i it is convenient to associate two index 
sets: 

Jnd(A {0) ) = {i G N : A^ > 0} (14a) 
Indi (A (0) ) = {i G N : Af } = -} (146) 

v p 

Let us present vector A^ G Wn, p ,i as the convex combination 

A(°) = P fi L » Wpll (L ) + (1 - p/i Lo )A (1) (15) 

where (see Eq.(10)) 

w P ii\ 



(ao = ^5>, ( 16 ) 



«gl ( , 



A (0) Lo A (0) 

A w = y A ; + y _\ (17) 



ieL ieN\L 
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and require the residual vector A^ to be representable. This requirement 
imposes the following restrictions on the admissible values of parameter jj, L °: 



o < ; ~1° <±, ieL 

— l—pn^O — p' KJ 



< <-, ie N\L Q 

The frontier solution of system (18) is 



(0) 



(18) 



/j, L ° = min{min{A- Uj }, min { 



1 



€N\L p 



(19) 



If n L ° ^ then we arrive at non-trivial representation of diagonal A^ as a 
convex combination of vertex w p i\(Lq) and a certain representable residual 
vector A«. From Eq.(19) it is easy to see that the additional condition 
fi L ° ^ holds true if and only if subset L satisfies the restriction 



Jnrfi(A (0) ) C L C Ind{\ (0) ) 

V 

Iterating of Eq.(15) leads to the following expression 



(20) 



fc-i 



i=Q 



.3=0 



PH Li w pil {Li) + 



'fc-i 



H(l ~P^) 



i=Q 



where 



and 



H u = min{min{A ; W } min - A ; w }} 

leLt ieN\Li p 



A (fe) (21) 

(22) 
(23) 



for % — 0, 1, . . . , k — 1. 

Definition 3. Sequence (L , Li, . . . , Lj, . . .) of p-element subsets of iV is 
called A-admissible if for each % = 0, 1, . . . subset Lj satisfies the condition 
(23). 

Theorem 2. For any vector A*- -* G Wjv, p ,i the residual vector in iteration 
formula (21) vanishes after a finite number of steps. 

Proof. First let us note that the number of nonzero components of 
representable residual vector A^ can not be less than p. If this number 
is equal to p then A^ just coincides with the vertex w P ii(Lk) where = 
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Ind(\( k >), and the residual vector A^ fe+1 ^ vanishes. Let us suppose that the 
number of nonzero components of \^ is greater than p. From Eqs.(15), 
(17), and (19) it readily follows that there exists index G Ind^X^) such 
that Af is necessarily equal either to zero or to -. To complete the proof 

it is sufficient to show that if A„-^ = - then A,- fc+1 ' ) = -. Condition (23) 
implies that all the indices i G N such that Af = - should belong to L k 
because in the opposite case the parameter \i Lk would be equal to zero. If 

^ = Af ) > and Af } = ± then Af +1) = i^fr = i. If, on the other hand, 



L 



=0 



Vi i Li e Li (24) 



f> = ±-Xg> > and Af > = I then 1-p/^ = P Af } and A? +1) = ^ = I 
□ P '* 

Corollary 1. The set of solutions of the Coleman's system (12) is the 
convex hull of ( n ) vertices w p ±i(L). 

Corollary 2. The number of vertices in expansion of a given density 
diagonal obtained on the base of the recurrence formula (21) is not greater 
than the number of its components different from zero. 

Corollary 3. A-admissible sequence (Lo, L\, . . . , L kx ) generated recur- 
rently on the base of the iteration formula (21) includes pairwise distinct 
p-element subsets and 

"i-l 

A(£o, Li, . . . , L kx ) = ^ 

8=0 

is a diagonal of p-electron density matrix such that 

i C P- 1 A(L ,L 1 ,...,L fc J = A. (25) 
pi 

It is to be noted that Theorem 2 is just a specification of the fundamental 
theorem by Caratheodory [TT] : 

Theorem 3. Let X C W 1 . Then any vector x G Conv(X) may be 
presented as a convex combination of no more than n + 1 vectors from X. 
Modern proof of this result may be found in • 

From Corollary 3 if follows that any mapping A — > s\ where s\ is a A- 
admissible sequence compatible with the iteration formula (21) determines 
some global section (right inverse) 7Tif P of the contraction operator ^c p ~ 1 
that is the mapping from Wn, p ,i to Tn }P such that 

j^-^piX) = A (26) 
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for any A G W^, p ,i- As it is seen from Eq.(22), sections constructed on the 
base of the recurrence relation (21) are not linear and even not differentiable 
in classic sense. 

The most ambitious task arising in the frameworks of the approach out- 
lined is to try to develop efficient methods for direct optimization of energy 
as a function of diagonal of the first order density matrix. General scheme 
embracing the whole class of such methods may be described as follows. 

1. Some section(s) of the contraction operator should be chosen. 

2. Using available section, it is possible of associate with some trial diag- 
onal A ensemble of p-electron determinant states and to determine squares 
of the CI coefficients: 



ICV 



|2 



'i-1 
j=0 



W Li (27) 



(see Eq.(24)). 

3. Construct average energy 



k 



A 



E x (fa,fa, • • • Ak x ) = J2 cos ^ - ^)|ClJ|C l .| < Li\H\Lj > (28) 

i,j=0 

as a function of phases fa. 
4. Minimize the function 

E^-.X^mmEM (29) 

to determine optimal diagonal and its expansion via vertices w P ii(L). 

There are no serious problems in implementation of steps 2-4 of this 
scheme and the only complicated step is reasonable selection of mapping(s) 
7Ti| p (note that in general several different sections may be employed in the 
course of the energy optimization). It is rather difficult to estimate a priori 
the quality of some chosen concrete section 7Ti| p . There are two readily 
coming to mind general algorithms to construct such sections. Both of them 
involve full sorting of p-electron subsets of the spin-orbital index set. 

1. Maximization of parameter (22) on each iteration: On the k-th step 
current may be determined from the condition 

a Lk = max <min{mm{Ap^}, min { ApH) 1 . (30) 

l { leL 1 ieN\L p 1 JJ J 
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In this case it is not necessary to take into account Eq.(23) explicitly. 

This section is probably optimal from formal mathematical viewpoint 
but has no physical idea behind it. Computer experiments show that in 
restoration process of such type high order excitations from the HF state 
contribute mostly. Even if exact FCI occupancies for the ground state are 
chosen, the restoration produces ensemble of determinant states that involves 
HF determinant and excited determinants that practically do not interact 
with the HF one. 

2. Energy minimization: On the k-th iteration among subsets satisfying 
the condition (23) it is chosen the subset such that the lowest eigenvalue 
of p-electron Hamiltonian in the basis {|L >, \L\ >,..., \L^ >} is minimal. 

This is undoubtedly the best possible section of the contraction operator. 
Unfortunately, the use of this section for the energy minimization is of no 
sense because it is equivalent to a certain CI scheme that can be described 
as follows. 

1. First it is necessary to fix the maximal number mdet of determinants 
in wave function expansion and put k = 0; 

2. Put k — k+1; 

3. Sort all determinants different from the already chosen and select the 
one that corresponds to the lowest eigenvalue of the Hamiltonian in the basis 
of k determinants {\L\ >, \L 2 >,..., >}. If k < mdet, return to step 2. 

Finally in a certain sense optimal basis involving not greater than mdet 
determinants will be obtained. This scheme is based on the well-known 
bracketing theorem of matrix algebra (see, e.g..|13j) and is used in quantum 
chemistry for years in different modifications to select initial determinant 
space for multi- reference CI calculations ■ in our opinion this scheme 

is interesting in its own right as a self-sufficient one when a relatively small 
number of leading determinants should be constructed from active orbitals 
with close orbital energies (the case that occurs extensively in transition 
metal complexes) because 

(1) CI spaces of huge dimensions can be efficiently handled and disk 
memory usage is minimal; 

(2) Calculations can be easily restarted; 

(3) Algorithms are trivially parallelized and if, say, PC clusters are used, 
data transfer via local net is minimal; 

(4) It is easy to handle both single excited state and a group of successive 
states. 

For the restoration purpose the above scheme can be considered as a 
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certain benchmark one because it gives the best possible occupancies and 
energy that can be obtained on the basis of the restoration routine described 
by Theorem 2. 



Conclusion 

General theorem establishing a connection between diagonal of the first 
order density matrix and a certain set of many-electron wave functions is 
proved. It is shown that rigorous energy expression involving only one- 
electron density becomes well-defined as soon as a certain right inverse of 
the contraction operator is chosen. For a fixed representable diagonal of the 
first order density matrix there exist quite a number of ways to restore p- 
electron determinant ensembles that are contracted to the diagonal under 
consideration. Each such way is in fact a path of a rather complicated graph 
with its vertices labelled by admissible (in sense of definition 3) p-element 
spin-orbital index sets. The main problem arising in implementation of op- 
timization schemes based on such energy expressions is the lack of general 
simple algorithms for selection of admissible paths for restoration of wave 
functions from one-electron densities. Such algorithms, besides requirement 
being simple, should generate paths close in a certain sense to ones obtained 
by the benchmark calculations based on the bracketing theorem. Search for 
such algorithms is in progress now. Note in conclusion that the recurrence 
formula (21) can be easily generalized to treat densities of higher order and 
the only obstacle here is the lack of the complete set of inequalities for ana- 
lytic description of the polyhedron W^,p, q in the case q > 1. 
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